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O . ABSTRACT 

We study afterglow flares of gamma-ray bursts (GRBs) in the framework of 
^ ! the late internal shock (LIS) model based on a careful description for the dynam- 

ics of a pair of shocks generated by a collision between two homogeneous shells,. 
First, by confronting the model with some fundamental observational features 
^ ! of X-ray flares, we find some constraints on the properties of the pre-collision 

shells that are directly determined by the central engine of GRBs. Second, high 
energy emission associated with X-ray flares, which arises from synchrotron self- 
^ ! Compton (SSC) emission of LISs, is investigated in a wide parameter space. The 

\Q ', predicted flux of high energy flares may reach as high as ~ lO^^erg cm~^s~^, 

^ ' which is likely to be detectable with the Large Area Telescope (LAT) aboard the 

Fermi Space Telescope (formerly GLAST). 

OO . 

Q , Subject headings: gamma rays: bursts — radiation mechanism: nonthermal 

> 

. 1. Introduction 

Since the launch of the Swift satellite in 2004, the X-ray telescope (XRT) aboard has 
revealed several new features of early X-ray afterglow light curves of gamma-ray bursts 
(GRBs) (e.g., Mezaros 2006; Zhang 2007 for reviews): (i) the transition from prompt phase 
to afterglow phase usually exhibits steep decline of X-ray flux, which is generally interpreted 
as the high-latitude tail of prompt emission (i.e., curvature effect; Fenimore et al. 1996; 
Kumar & Panaitescu 2000;), (ii) a shallow decay usually follows the steep segment during 
the first hours, which could be caused by a spread of Lorentz factors in GRB ejecta (Rees 
& Msezaros 1998) or a continuous energy injection into GRB ejeta due to lasting activities 
of GRB central engines (e.g., Dai & Lu 1998a, 1998b; Zhang & Meszaros 2001a; Dai 2004; 
Fan & Xu 2006; Yu & Dai 2007), and (iii) bright X-ray flares superimposing on underlying 
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afterglow emission were observed from nearly a half of the Swift GRBs (Burrows et al. 2005; 
see Chincarini et al. 2007 and Falcone et al. 2007 for more recent statistic studies). These 
discoveries are of great importance for revealing the nature of GRB central engines, and 
especially, X-ray flares are widely accepted to be due to some delayed intermittent activities 
of the central engines. 

In the sight of model, the delayed intermittent activities might be episodic accretion onto 
a central object due to a chopped accretion disk (Perna ct al. 2006), or episodic accretion 
due to a modulation of the accretion flow by a magnetic barrier (Proga & Zhang 2006), 
or magnetic rcconncction on a nascent differentially-rotating massive neutron star (Dai et 
al. 2006), etc. Many scenarios were designed in the past few years, but none of them is 
conclusive (see Zhang 2007 for a review). However, the key point is that, in almost all of 
these scenarios, internal dissipations in late-ejected materials are usually required to produce 
flare emission. The leading mechanism for internal dissipation is late internal shocks (LISs; 
Burrows et al. 2005; Fan & Wei 2005), which are generated by collisions among some shells 
ejected at different times with different Lorentz factors and energies. Hitherto, the properties 
of LlS-produccd emission have been discussed for a few times (e.g.. Fan & Wei 2005; Wu 
et al. 2005; Zhang et al. 2006; Zou et al. 2006; Fan et al. 2007; Galli & Guetta 2008). 
However, the properties of the pre-coUision shells, which are directly determined by the late 
activities of the central engines, were not concerned in their works. In this paper, in contrast, 
we give a more detailed description for the dynamics of LISs including internal forward and 
reverse shocks, and then investigate the predicted X-ray flare luminosity and the shape of 
light curves within a wide model parameter space. By confronting the theoretical predictions 
with some observational features of X-ray flares, we can find some constraints on the model 
parameters. 

Furthermore, it must be helpful to constrain the model more stringently by performing 
a simultaneous observation in the high-energy 7-ray bands for X-ray flares in the Fermi 
era. As suggested previously by Wang, Li, & Meszaros (2006), Fan & Piran (2006), and 
Galli & Guetta (2008), high energy flares are expected to arise from inverse Compton (IC) 
scattering of low-energy flare photons off some relativistic electrons. These electrons belong 
to either GRB ejecta (i.e. external inverse Compton, EIC) or flare ejecta self (i.e. synchrotron 
self-Compton, SSC). Based on the same consideration, we discuss the detectabihty of high 
energy flares in the Fermi era with the careful description for LIS dynamics along with the 
constraints from X-ray observations. 

This paper is organized as follows: in section 2, we describe the dynamics of a pair of 
shocks arising from a collision between two homogeneous shells and the resulting synchrotron 
radiation of the shocked electrons. In section 3, the model is tested by some observational 
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features of X-ray flares. In section 4, we estimate the accompanying SSC high-energy emis- 
sion and its detectabihty in the Fermi era. Finally, a summary is given in section 5. 



2. The model 

After the prompt phase of a GRB, the central engine might still be able to eject some 
separate shells with different Lorentz factors and energies due to its delayed intermittent 
activities. CoUisions between these shells generate LISs, which give rise to the observed 
X-ray flares. For simplicity, the scope of the paper is restricted within one collision between 
two homogeneous shells. 



2.1. Dynamics 

We set the time zero point at the GRB trigger and measure time in the observer's frame. 
At a time of iej,A, the central engine ejects a shell denoted by A, which moves at a constant 
bulk Lorentz factor 7a and carries an isotropic-cquivalent kinetic energy luminosity Liso,A- 
The particle number density of the shell measured in its comoving frame (denoted by a prime 
hereafter) can be calculated by 

/ -^iso,A /-,\ 

where R is the radius of the shell. Some time (Atgj) later, another shell denoted by B with 7b 
and Liso.B is assumed to be ejected again. We require 7b > 7a in order to let shell B catch up 
and coUide with the prior shell A. Consequently, a coUision between A and B takes place at 
the radius i?coi = PAPBcAtej/iPB - Pa) and the time ^coi = ^ej,A + /3AAiej/(/3B - /3a) - -Rcoi/c 
For (7a, 7b) > 1, they read 

^ 27icAt, 

^°'"1-(7a/7b)^' ^ ^ 

^col — ^ej.A ~l~ Z 7 7 — ^onset- (3) 

1 - (7a/7b)2 

The collision time tcoi determines the observed onset time tonsct of a fiarc in physics. Strictly, 
tonset should be mildly larger than tcoi since the physical onset is usually buried by normal 
afterglow emission. 

Due to the collision, a pair of shocks are generated, including a forward shock and a 
reverse shock that propagate into shells A and B, respectively. Separated by the two shocks 
and a contact discontinuity surface, the system is divided into four regions: (1) unshocked 
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shell A, (2) shocked shell A, (3) shocked shell B, and (4) unshocked shell B, bulk Lorentz 
factors of which are 71 = 7a, 72 = 73 = 7, and 74 = 7b. Considering the motion of the 
shocked regions relative to unshocked regions 1 and 4, respectively, two relative Lorentz 
factors can be calculated by 

1 /71 , 7 \ 1 ^7 , 74\ 

721 = - — + — , 734 = - — + — • (4) 

2 V7 71/ 2 V74 7/ 

Then, according to the jump conditions between the two sides of a shock (Blandford & 
McKee 1976), we can calculate the internal energy densities of the two shocked regions by 
62 = (721 — l)(472i + 3)n[mpC^ and Cg = (734 — 1)(4734 + 3)n'^'mpC^, where n[ = n'^ and 
n'^ = n'g. The mechanical equilibrium between the two shocked regions requires e'2 = 63, 
which yields 

(72i-l)(472i + 3) _< _ fLA 

(734 -l)(4734 + 3) n[ \lJ\^J ^' 

For four given parameters, Li, L4, 71, and 74, that describe the basic properties of the pre- 
collision shells, we can obtain the value of 7 by solving equations (jlj) and In particular, 
for four limits shown in the (74/71, L4/L1) parameter space (see Figure 1), these equations 
can be solved analytically. For 74 ^ 71, we have (i) if L4/L1 ^ (1/7) (74/71)"^, 

2 

721 = 7?^ » 1 & 734 - 1 ^ 3^ = e « 1; 7 = 74(1 - V2e), (6) 
271 'III 

which means the forward shock is relativistic and the reverse shock is Newtonian (RFS & 
NRS); (ii) if 16 « L4/L1 « (1/16) (74/71)', 

rl/i 1/2 1/2 

/ 74 1 B 74 1 /1/4 1/2 1/2 

721 = f75- > 1 & 734 = ^^^^ 1/2 7 = / ^ 7/ 74 , 7 

which means both the forward and reverse shocks are relativistic (RFS & RRS); (iii) if 
LJL, < 7, 

721 - 1 ^ (4 = ^ « 1 734 = » 1; 7 = 7i(l + V^), (8) 
771 271 

which means the forward shock is Newtonian whereas the reverse shock is relativistic (NFS 
& RRS). Finally, (iv) for 74 ^ 71, both the forward and reverse shocks are Newtonian 
(NFS & NRS). Since 71, 74, and the ratio / are unchanged with the moving of the shells, 
the values of 7 as well as 721 and 734 can keep constant before the shocks cross the shells. 
In principle, when one shock crosses the corresponding shell, the other shock should be 
decelerated in the rest frame of its upstream material. However, if the energy carried by 
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the remaining unshocked material is much less than the energy of the total shocked material 
at that time, 7 would be not changed significantly. So, we can treat with 7 as a constant 
always approximately if we assume the crossing times of the two shocks to be equal more or 
less. After both shocks vanish, the merged shell no longer interacts with any other materials 
until it is caught up with by the third late-ejected shell. 

The above analysis on kinematics and dynamics yields the evolution of the radius of the 
system after the collision. 



where we redefine the time T by resetting the time zero point from the GRB trigger to the 
fiare onset time. We can define an initial expansion time from equation ([9]), 



and find that the increase of the radius before Texp can be ignored (i.e. R ~ constant) but 
after Tgxp the radius increases linearly with time (i.e. R gcT). Moreover, as the propagation 
of the shocks, the total electron numbers can be calculated by A^e,2 = ^T^R^TT''i{.l2i(32ihP)l'^cT 
and Ne^z = ^''^R'^n'^{l34:l334:/lP)'l'^cT for the two shocked regions (Dai & Lu 2002), both of 
which are proportional to T before the shock crossing. 



For the two shocked regions, fractions €3 and eg of the internal energy are assumed to be 
carried by magnetic fields and hot electrons, respectively. Then, the strength of the magnetic 
fields reads B'- = (STre^e^)^/^, whose variation is determined by the evolution of e[. Denoting 
the comoving volume of a shocked region as V/, we adopt oc V'^'^ oc R~^ by ignoring a 
possible spreading of the hot shocked materials. In the presence of the shocks that might 
be able to suppress the spreading of the hot materials, the neglect of the spreading may be 
plausible. However, after the shock crossing time T^s (we assume here the two shocks cross 
at a similar time), the spreading of the hot materials into the vacuum cannot be ignored 
and the hot materials should experience an adiabatic cooling. During this phase, we assume 
that the volume of the merged shell is determined by a simple power-law as V/ oc R^, where 
s is a free parameter and its value is taken from 2 to 3. Then, the particle number densities 
would decrease as n- oc V'^ oc R~^ and the internal energy densities as e- oc 
due to adiabatic cooling. Therefore, a multi-power temporal behavior for the magnetic field 



R{t) = i?col + 2j^c{t - teol) = Real + 2j^cT, 



(9) 




(10) 



2.2. Synchrotron radiation 
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strength can be found: 



, T°, T < Tex ; 

T -'- ^ -'- exp ) 



TO, T < Texp; 
oc T-\ Texp < T < T„s; T,,, > Te^p. (12) 

T^— 2s/3 T ^ T ■ 
-'- 1 -'- -'- crs ) 

As usual, we assume that electrons are accelerated by the shocks to distribute as 
dngj/d7'g j oc for 7'g j > 7'e,m,j' where the minimum Lorentz factor is defined as 

l'e,m,i = eeCp{mp/me) (Trei"!) with Cp = {p-2)/ (p-l) and 7rei = 721 or 734. Additionally, the 
cooling Lorentz factor, 7'ecj = 67rmec/[(l +yj)crT-B'^7T], also needs to be introduced, above 
which the electrons lose most of their enegies. The Compton parameter Yi defined as the 
ratio of the IC to synchrotron luminosities can be estimated by ^ [{Arfiee/eB + 1)^^^ — l]/2 
with rii = min[l, (7e,c,i/7e,m,i)^~^] (Sari & Esin 2001). Then, by calculating two characteristic 
frequencies and a peak flux density, 

"-=2^''^^^'-^^' "^^--'^4^^^^^^' ^'^^ 

where di is the luminosity distance of the source, a multi-power synchrotron spectrum con- 
tributed by a shock can be written as (Sari et al. 1998) 



1/3 
(</-l)/2 



(^) , n<i^<i^h-, (14) 



(q-l)/2 ^ ^ -p/2 



where vi = min(i/„_j, i/c,i), = max(i/^,j, i/^j), and g = 2 for v^^i < Vrn,i and g = p for 

To obtain X-ray light curves further, the temporal dependence of the characteristic 
quantities also needs to be presented as follows: 



(i) In the case of T^rs < Texp, 



TO T < T ■ 

) 1 ^ -'expi 

^m,i OC \ rp-2s/3 T ^ T ^ ^ 

-'- 1 -'- -'- exp ) 

^c,i OC <i ^2.-2^ ^ ^ ^r^2l ^^^^ 
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1 J. \ J. crs ) 

Fi/^xaax.,i ^ \ ^ i ^crs < < Texpj (-'-''^) 

r-2V3^ r>rexp; 



(ii) In the case of Tcrs > Texp, 



T < Texp; 

i^m,iOC<{ T-l, rexp<r<rers; (18) 
7^— 2s/3 T ^ T 

1 ^ ^ ^ crs ) 

T <T ■ 

) J. \ -I. exp ) 

i/c,i OC <( T, Texp < T < Tcrs; (19) 

y2.-2^ r>Ters, 



T < T ■ 

J- ^ -I- exp ) 




-^i/,max,j \ ) -^xp T <i Tf-j-g, (20) 



From the above expressions, we can see that i/c,i reaches its minimum value at Tcxp, while 
Urn^i starts to decrease at the same time. So, the relationships between Vc^i and Vrn,i as well 
as the spectra and light curves can be found easily by fixing and comparing the values of v^^i 

and Vrn,i at Texp. 



3. X-ray flares 

Observations have shown that X-ray flares may consist of one or a few pulses. We 
consider one pulse to be mainly produced by one collision between two shells. Here we do 
not try to fit X-ray observational data in detail, which requires a more elaborate model that 
takes the fine structure of pre-coUision shells into account. Instead, we only test the model 
by some fundamental observational features of X-ray fiares. 



3.1. X-ray luminosity 

Statistic studies found that the average flare fluence (in 0.2 — 10 keV band) is approxi- 
mately a factor of ten less than the fluence of prompt emission as ~ 10~^ erg cm~^ (Falcone 
et al. 2007) and the flare peak times t-p^ak concentrate into the range from 100 to 1000 s 
(Chincarini et al. 2007). Moreover, according to St/tpeak ^ 1, the temporal width of flares. 
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6t, can be estimated to be from a few to several ten seconds. By considering the spread of 
the distributions of these quantities, we suggest to take ~ 10~^ erg s~^ cm~^ as the lower 
limit for the peak flux of X-ray flares. For a typical luminosity distance = 10^^ cm for 
GRBs, we give ~ 10^^ erg s~^ as the lower limit for X-ray flare luminosity. 

In order to derive some constraints on the model from this luminosity lower limit, first we 
show two (71, 74) parameter spaces in Figure 2 with certain values of Li and L4: L4/L1 ^ 1 
for the left panel and L4/L1 = 1 for the right panel, which correspond to different dynamic 
cases. According to the variation of and i/* jl| the parameter spaces can be roughly divided 
into four regions denoted by 'a, b, c, d', where different relationships among z/j^, z/*, and 
i'x{= 2.4 X 10^''Hz) are given as listed in Table 1. We then use some shaded contours to show 
the regions where model-predicted X-ray luminosity exceeds the observational lower limit. 
Obviously, region 'a' is restricted significantly by this luminosity constraint due to the slow 
cooling of electrons (as indicated by z/j^ < z/x < t'*) that leads to a low radiation efficiency. 
To be more specific, we can conclude further that (i) in order to produce sufficiently strong 
flare emission, it is required that 74 > few x 71. This means that at least one of the 
forward and reverse shocks is mild-relativistic; and (ii) when 74 ^ few x 7^, we can obtain 
the highest theoretical X-ray luminosity, because this parameter region roughly locates the 
boundary between regions 'c' and 'd', where the relationship u* < ux ^ can be found, 
indicating a high X-ray radiation efficiency. Second, in order to uncover the dependence of 
the emission on parameters Li and L4, we also show two corresponding parameter spaces 
in Figure 3 with certain values of 71 and 74. A minimum value of ~ 10^° — 10^^ erg s~^ for 
L4 can be found, but the constraint on Li is loose. This indicates that the resulting X-ray 
luminosity is mainly determined by the kinetic-energy luminosity of the lagged rapid shell 
rather than that of the leading slow shell. 

Finally, in the above calculations for X-ray luminosity, we ignored possible synchrotron 
self-absorption of the X-ray photons, which is able to suppress the X-ray emission. The 
synchrotron self-absorption thickness at the X-ray band can be calculated by (Panaitescu & 
Kumar 2000) 



^The star superscription represents that the values are calculated at the time of Tcxp- We obtain the 
values of and ly* from the spectrum due to the shock that dominates the X-ray emission and thus the 
subscription i of the quantities is omitted in this section. 




(21) 
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where ui = min(i/^, i/c), Uh = max(i/^, i/c), 7^^^ = mm(7g 7^^), and g = 2 for i/c < i^m and 
q = p for i/f. > I'm- By scanning the whole parameter spaces shown in Figures 2 and 3, wc find 
that the values of Tssa,vx always much lower than unity due to the large internal shock 
radius, which means the synchrotron self-absorption of the X-ray photons can be ignored 
safely. 



3.2. Shape of X-ray light curves 

We here test the model by considering the shape of the observed X-ray flare light curves, 
a basic feature of which is the rapid rise and fall. Following the frequency relationships given 
in Table 1, we can find three possible types of the theoretical X-ray light curves, as shown 
schematically by the black lines in Figure 4 for Tcrs — Tcxp where the time zero is set at the 
flare onset time. The light curves break at several characteristic times and the corresponding 
temporal indices {a — —d\ogFy/d\ogT) are listed in Table 1. All of the segments after Tcrs 
are steepened significantly by the shock crossing effect. By confronting these theoretical 
light curves with the observed ones, on one hand, it is easy to understand the observed 
steep rise by resetting the time zero point of the log-log figure from the fiare onset time 
to the GRB trigger. On the other hand, as discovered by Liang et al. (2006), the rapid 
decline of most X-ray fiares seems to be consistent with the curvature effect that predicts 
a temporal index being equal to the simultaneous spectral index plus 2. As shown by the 
black lines in Figure 4, the intrinsic decline slope of the last segment is a = {sp H- s)/3 for 
all types of the theoretical light curves and the corresponding spectral index is (p — l)/2 due 
to < < ^c- For s — 2) and 2 < p < 3, the inequation 

^>V^ + 2 (22) 

can be satisfied easily. This inequation indicates that, in any case in our model, the observed 
X-ray fiux decay should be dominated by the curvature effect, as shown by the grey lines 
in Figure 4. In addition, although some relatively flat segments with a — (3 — s)/3 or 
{sp — 2s -I- 3)/3 appear in types II and III light curves, we still cannot rule out any one of 
them absolutely, because these segments that could not be very far from the flare onset time 
may be steepened by the time zero effect dramatically. 



4. High energy flares 



Because a part of the synchrotron photons would be boosted to higher energy by their 
IC scattering off some relativistic electrons, high energy counterparts of X-ray flares are 
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expected naturally. In the following calculations, we restrict our attention within the LIS 
model and SSC scenario. Before a specific calculation for the SSC emission, we first estimate 
the flux sensitivity of the LAT aboard the Fermi by 



thr 



1.33 X 10-9 ( (-^) erg cm^^ s-\ (23) 



following Zhang & Meszaros (2001b) who adopted the criterion that at least five photons 
are collected if the instrument is source dominated, where E is the photon energy and t 
is the integration time. The effective area A^g of the instrument is taken as a constant of 
6000 cm^ and the dependence of the area on the photon energy is neglected. However, after a 
transition time of ~ 2.4 x lO^s (Gou & Meszaros 2007; Yu et al. 2007), the sensitivity should 
start to scale as due to a limitation by the background for a long-time observation. 

Following Sari & Esin (2001), the SSC spectrum contributed by a shock can be obtained 
by shifting the seed synchrotron spectrum to higher energy range, i.e. estimating the two 
break frequencies and the peak flux density of the SSC spectrum respectively by 



z/^. = 2V' z/ • z/^ = 2V' u ■ . = ^lI^F (24) 



This approximative treatment is finely valid when the effect of the Klein-Nishina suppression 
is unimportant. Considering the highest-energy electrons whose energy enter the Klein- 
Nishina regime, we refer to the third break in the SSC spectrum, 

2 2 4 

max(z/„_i, z/^j) 

above which the SSC spectrum follows F^^i oc u~^p~^^ (Gupta & Zhang 2007; Fragile et al. 
2004). Therefore, the approximative SSC spectrum can be summarized as follows (Gupta & 
Zhang 2007) 



, V3 

u 



F'^ ■ = F'^ . X < 



-(q-l)/2 / N -p/2 



(26) 



Vh<V < Z/KN,i, 
Ul J \ I^h J \i^KN ' ' ' 

where i^i = min(z/^ i/^), Uh = max(i/^_^, z/^), and g = 2 for z/^ < z/^^^ and g = p for 
z/^ ■ > z/^ .. 

c,t m,i 



Using equation (|26|) . we calculate the peak flux of GeV 7-ray flares as shown by the 
shaded contours in Figure 5. Comparing these GeV emission contours with the dashed 
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contours that correspond to the X-ray luminosity, a positive correlation between these two 
emission components can be found. For the GRBs at a typical distance of lO^^cm, the 
high energy counterparts of the relatively brighter X-ray flares could be detected by the 
LAT, whereas those associated with weaker X-ray flares leak. The LAT sensitivity here 
(few X 10^^ erg s~^ cm~^) is calculated for a typical flare onset time of several hundred 
seconds and represented by the dash-dotted line in the figure. Furthermore, according to 
the expression of ?a [(Arjiee/eB + 1)^''^ — l]/2, we know that the relative importance of the 
SSC and synchrotron emissions may be sensitive to the parameters and eg. Therefore, we 
show the GeV 7-ray flux at the time of Texp varying in the (eg, €3) parameter space in Figure 
6. It can be seen that the high energy flux is mainly sensitive to and equipartition values 
for Ce are required. In addition, we would like to show some example spectra numerically in 
Figure 7 using a more elaborate code that was developed in Yu et al. (2007). 



5. Summary and discussion 

The LIS model is usually employed to explain the observed GRB afterglow X-ray flares. 
However, a careful description for LIS dynamics and some observational constraints on it 
still need to be investigated. Based on this consideration, we studied the properties of LIS- 
produced emission in the framework of a simplified paradigm, i.e. internal forward-reverse 
shocks generated by a collision between two homogeneous shells. With a lower limit for 
the observed X-ray luminosity, we found a mildly high ratio of the Lorentz factors of pre- 
coUision shells, which leads to at least one mildly-relativistic internal shock. Our results 
also show that the brightest X-ray flares might imply a high variability of Lorentz factors as 
indicated by 74 = few x jf. The luminosity of the X-ray flares is mainly determined by the 
kinetic energy luminosity of the delayed rapid shell rather than the leading slow shell. After 
an investigation of the characteristic frequencies in a wide parameter space, three types 
of theoretical X-ray light curves arc found, all of which are ended by a very steep decay 
with a ^ p + 1. This indicates that the flare emission during the decay phase is probably 
dominated by the curvature effect, which is consistent with the observational inference found 
by Liang et al. (2006). 

We also investigated the peak flux of the GeV 7-ray counterparts in the SSC scenario. 
By comparing the GeV flux with the flux sensitivity of the Fermi LAT, we found that the 
high energy flares associated with relatively brighter X-ray flares could be detected by the 
LAT for a distance of lO^^cm to the source, where an equipartition value of is required. 
This possible detection will be very helpful to discriminate different origins of high energy 
flares and different models for X-ray flares. As mentioned above, two types of high energy 
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flares are predicted by the LIS model including the ones due to the SSC and EIC emission. 

In the SSC case, a good temporal correlation between the X-ray and high energy flares are 
expected, whereas a signiflcant temporal extension appears for high energy flares in the EIC 
case (Fan et al. 2007). In addition, besides the LIS model, some authors suggested that 
X-ray flares may be produced by a delayed external shock and the corresponding high energy 
emission was also expected (Galli & Piro 2007). But for the delayed afterglow model it is 
difficult to explain the reoccurrence of X-ray flares in one GRB afterglow. 

Acknowledgements 

This work is supported by the National Natural Science Foundation of China (grants no. 
10221001, 10640420144, and 10873009) and the National Basic Research Program of China 
(973 program) No. 2007CB8 15404. YWY is also supported by the Scientiflc Innovation 
Foundation of Huazhong Normal University, the Visiting PhD Candidate Foundation of 
Nanjing University and the National Natural Science Foundation of China (grant 10773004). 

REFERENCES 

Blandford, R. D., & McKee, C. F. 1976, Phys. Fluids, 19, 1130 

Burrows, D. N., Romano, P., Falcone, A., Kobayashi, S., Zhang, B., Moretti, A., O'Brien, 
P. T., Goad, M. R. et al. 2005, Science, 309, 1833 

Chincarini, G., Morettil, A., Romano, P., Falcone, A.D., Morris, D., Racusin, J., Campana, 
S. et al. 2007, ApJ, 671, 1903 

Dai, Z. G. 2004, ApJ, 606, 1000 

Dai, Z. G., & Lu, T. 1998a, A&A, 333, L87 

Dai, Z. G., & Lu, T. 1998b, Phys. Rev. Lett., 81, 4301 

Dai, Z. G., & Lu, T. 2002, ApJ, 565, L87 

Dai, Z. G., Wang, X. Y., Wu, X. F., & Zhang, B. 2006, Science, 311, 1127 

Falcone, A. D., Morris, D., Racusin, J., Chincarini, G., Moretti, A., Romano, P., Burrows, 
D. N., Pagani, C. et al. 2007 ApJ, 671, 1921 

Fan, Y. Z., Piran, T., Narayan, R., & Wei, D. M. 2008, MNRAS, 384, 1483 



Fan, Y. Z., & Piran, T. 2006, MNRAS, 370, L24 
Fan, Y. Z., & Wei, D. M. 2005, MNRAS, 364, L42 
Fan, Y. Z., & Xu D. 2006, MNRAS, 372, L19 

Fenimore, E. E., Madras, C. D., & Nayakshin, S. 1996, ApJ, 473, 998 

Fragile, R C, Mathews, G. J., Poirier, J., & Totani, T. 2004, Astropart. Rhys., 20, 591 

Kumar, P. & Panaitescu, A. 2000, ApJ, 541, L51 

Liang, E. W. et al., 2006, ApJ, 646, 351 

Gain, A., & Guetta, D. 2008, A&A, 480, 5 

Gain, A., & Piro, L. 2007, A&A, 475, 421 

Gou, L. J., & Meszaros, P. 2007, ApJ, 668, 392 

Gupta, N., & Zhang, B. 2007, MNRAS, 380, 78 

Mezaros, P. 2006, Rept.Prog.Phys., 69, 2259 

Panaitescu, A., & Kumar, P. 2000, ApJ, 543, 66 

Perna, R., Armitage, P. J., & Zhang, B. 2006, ApJ, 636, L29 

Proga, D., & Zhang, B. 2005, MNRAS, 370, L61 

Rees, M. J., & Msezaros, P., 1998, ApJ, 496, LI 

Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 

Sari, R., & Esin, A. A. 2001, ApJ, 548, 787 

Wang, X. Y., Li, Z., & Meszaros, P. 2006, ApJ, 641, L89 



Wu, X. F., Dai, Z. G., Huang, Y. F., Feng, L. L., & Lu T. 2005, arXiv: |astro-ph/05 12555 

Yu, Y. W., & Dai, Z. G. 2007, A&A, 470, 119 

Yu, Y. W., Liu, X. W., & Dai, Z. G. 2007, ApJ, 671, 637 

Zhang, B. 2007, ChJAA, 7, 1 



- 14- 

Zhang, B., Fan, Y. Z., Dyks, J., Kobayashi, S., Meszaros, P., Burrows, D. N., Nousek, J. A., 
& Gehrels, N. 2006, 642, 354 

Zhang, B., & Meszaros, P. 2001a, ApJ, 552, L35 

Zhang, B., & Meszaros, P. 2001b, ApJ, 559, 110 

Zou, Y. C, Dai, Z. G., & Xu, D. 2006, ApJ, 646, 1098 



This preprint was prepared with the A AS IM^jX macros v5.2. 

Table 1: Four possible relationships (a, b, c and d) between u^, u*, and Ux and three types 
(I— III) of the corresponding X-ray flare light curves (characterized by the temporal indices). 



Regimes 


Types 


Temporal Indices (a) 


a: 


iy^<iyx< K 


I 


1 sp+3 
^' 3 


b: 


< < K < 


II 


1 fi sp— 2s+3 sp+3 


c: 




II 


1 n sp— 2s+3 sp+3 
i,U, 3 , 3 


d: 


ly* <iyx<iy*m 


III 


1 5 rv 3-s sp-2s+3 sp+3 
' 3' ' 3 ' 3 '3 
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Fig. 1. — Regions in the (74/71, L4/L1) parameter space where four hmits of the LIS dy- 
namics are given. 




Fig. 2. — X-ray luminosity in the (71, 74) parameter spaces. The regions where the model- 
predicted X-ray luminosity exceeds the observational lower limit (~ 10^^ erg s~^) are pre- 
sented by shaded contours, while the unshaded region is ruled out by this luminosity con- 
straint. Separating by solid lines, the parameter spaces are divided into different regions 
denoted by "a, b, c, d", where different relationships between u*, u^, and ux are given as 
listed in Table 1. The fixed values of the kinetic-energy luminosities of the shells satisfy 
L4/L1 ^ 1 for the left panel and L4/L1 = 1 for the right panel, and their corresponding 
dynamic cases can be found in Figure 1. The black region is forbidden due to 74 < 71. The 
other model parameters eg, eg, p, and A^ej are taken to be typical values of 0.03, 0.3, 2.5, 
and 100s, respectively. 
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Fig. 3. — X-ray luminosity in the {Li, L4) parameter spaces. The left and right panels cor- 
respond to relatively higher and lower values of the ratio 74/71, respectively. The meanings 
of the regions and the other model parameters are the same as those in Figure 2. 




Fig. 4. — Schematic illustration of theoretical X-ray flare light curves for Tcrs = Tcxp (black 
lines; the temporal indices of all segments are listed in Table 1). The curvature effect is 
exhibited by the grey hnes. The vertical dotted lines represent the break times of the hght 
curves, specifically, the time Tcm at which i/c — 1/^, the time (Tc) at which the break 
frequency 1/^ (i^c) passes through the X-ray band, the time Texp from which the radius 
increases linearly. 



Io9io^i 



Fig. 5. — The peak flux at T = Tgxp of GeV 7-ray flares due to SSC emission for GRBs at 
tlie luminosity distance lO^^cm in the (71, 74) parameter space. The regions where the GeV 
flux exceeds the Fermi LAT sensitivity (dash-dotted hne) are shown by shaded contours, 
while the SSC emission calculated in the unshaded region could be not detected by the LAT. 
The model parameters here are the same as those adopted in the left panel of Figure 2. To 
compare with the associated X-ray component, the X-ray luminosities are also shown by the 
dashed contours and labeled by logiQ[{uL^)x/ey^g s"^]. 
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Fig. 6. — Variation of the peak flux at T = Tgxp of GeV 7-ray flares due to SSC emission for 
GRBs at the luminosity distance 10^^ cm in the (ee, e^) parameter space. The contours are 
labeled by the values of logio[(i^-^i/)Gev/erg s~^cm~^]. The region where the GeV flux exceeds 
the Fermi LAT sensitivity (dash-dotted line) is shaded. The black region is forbidden due to 
es+Q > 1- The other model parameters are taken to be Li — 10^° erg s~^, L4 = 10^^ erg s""-*^, 
7i = 10, 74 = 300, p = 2.5, Aiej = 100s and iej.A = 400s. 
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Fig. 7. — An example numerically-calculated spectrum by combining the synchrotron and 
SSC spectra contributed by the two shocks for the luminosity distance of lO^^cm. The 
dash-dotted line denotes the Fermi LAT sensitivity. The model parameters are taken to be 
eB = 0.03, e, = 0.3, Li = 10^° erg s-^ L4 = lO^^ gj.g g-i^ ^ ^q, 74 = 300, p = 2.5, 
Atej = 100s and tej,A = 400s. 



